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摘要 在 基于 目标 轨道 运动 特征 的 预警 信息 推算 中 , 引入 高 程 信息 库 的 支持 , 使 对 地 球 的 假设 从 参考 椭 球 体 变 
为 更 真实 的 形态 ， 并 探讨 在 此 情况 下 提高 发 、 落 点 位 置 和 时 刻 的 预报 精度 的 方法 . 采用 数字 高 程 模型 进行 的 数 
值 实验 结果 表明 , 在 原 方法 的 基础 上 增加 少量 高 程 欠 代 计算 的 情况 下 即 可 提高 预报 精度 , 发 、 落 点 位 置 高 程 越 
- 高 的 地 区 精度 提高 越 明显 . 

一 关键 词 弹道 , 方法 : 数值 , 高 度 
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1 引言 fle，HRRP) 具 有 成 像 速 度 快 、 分 辩 率 高 的 特点 ， 
弹道 导弹 预警 在 航天 技术 和 戌 略 预警 中 具有 且 HRRP 能 提供 沿 雷达 视线 方向 强 散 射 中 心 的 分 
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课题 弹道 导弹 的 整 运行 过 程 分 为 助 推 段 、 中 段 和 ”目标 特征 库 的 模板 对 比 , 即 可 完成 目标 号 份 鉴别 , 
再 入 段 三 部 分 . 其 中 中 段 又 叫 自由 飞行 段 , 是 弹道 而 关于 发 、 落 点 信息 推算 的 研究 都 是 在 已 知 目标 
导弹 关闭 发 动机 后 在 大 气 层 外 飞行 的 过 程 , 这 个 阶 “为 弹道 导弹 的 前 提 下 , 利用 弹道 与 地 球 有 交点 这 一 
段 飞 行 时 间 最 长 , 占 整个 弹道 飞行 时 间 的 80% 以 上 ， ”特点 进行 推算 . 推算 方法 多 是 如 文献 由 介绍 的 使 用 
> 仅 受 自然 力 影响 , 其 轨道 特征 容易 被 掌握 . 因此 卡尔 曼 滤 波 或 改进 的 卡尔 曼 滤 波 对 目标 运动 方程 
C 中段 是 各 种 无 线 测量 设备 进行 捕获 、 跟 踪 、 身 份 。 ”进行 外 推 的 方法 进行 计算 . 这 类 外 推 方法 本 质 上 是 
c o0 ”和 鉴别、 计算 预 警 信息 (如 发 、 落 点 时 间 、 位 置 ) 的 主 。 “一 种 试探 法 , 即 每 一 次 外 推 结 束 后 需要 判断 是 否 ) 
© 要 阶段 . 入 地 球 内 部 , 如 果 进 入 地 球 内 部 则 需 采 用 二 分 法 反 
在 此 前 的 研究 中 , 弹道 导弹 中 段 的 身份 鉴别 和 复 外 推 最 终 求 得 结果 . 其 计算 效率 与 滤波 步 长 的 选 
发 、 落 点 信息 推算 是 两 个 不 同 的 课题 . 文献 [1-3] 取 密 切 相 关 , 最 合适 的 滤波 时 间 步 长 又 与 具体 的 弹 
等 都 基于 弹道 导弹 弹头 在 雷达 回 波 下 反映 出 的 物 A, 导致 很 难 有 普 适 的 最 优 步 长 可 供 选 择 . 此 
里 特征 ,例如 弹头 的 微 动 (如 进 动 、 章 动 、 雷 达 散 外 这 种 方法 未 考虑 目标 在 飞行 时 受到 的 摄 动力 , 也 
射 截面 (Radar Cross Section, RCS) 的 变化 ) 存 在 周 没有 加 入 高 程 信息 的 支持 . 也 有 文献 提出 把 导弹 轨 
期 性 特点 ， 可 通过 雷达 回 波 呈现 的 微 多 普 勒 效应 道 简化 为 二 体 模 型 下 的 椭圆 轨道 后 通过 几何 关系 
进行 分 析 和 识别 ， 又 例如 雷达 对 飞行 目标 进行 的 求解 弹道 与 地 球 交 点 是 , 但 这 种 方法 同样 没有 考虑 
高 分 辩 率 一 维 成 像 (High Resolution Range Pro- 摄 动因 素 和 高 程 的 影响 . 文献 [6] 则 另辟蹊径 , 根据 
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导弹 轨道 与 地 球 有 交点 , 而 卫星 轨道 与 地 球 无 交点 
这 一 特征 , 提出 了 一 种 同时 进行 快速 鉴别 目标 身份 
和 生成 预警 信息 的 方法 . 该 方法 计算 速度 快 、 处 理 
过 程 简 单 稳健 、 计 算 效 率 高 , 且 考 虑 了 地 球 对 目标 
的 主要 摄 动 项 . 但 在 使 用 该 方法 生成 预警 信息 时 ， 
地 球 被 粗略 地 简化 为 一 个 参考 椭 球 体 , 高 程 信息 依 
然 被 略 去 . 这 样 虽然 减少 了 计算 量 、 提 高 了 计算 速 
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(1) 


其 中 了 为 目标 于 t 时 刻 的 真 近 点 角 . 求解 (1) 式 的 方 
法 就 是 通过 轨道 运动 特点 鉴别 目标 身份 和 推算 预 
和 警 信息 的 核心 内 容 ，f 有 解 就 证 明 目 标 与 地 球 有 交 
点 , 通过 轨道 根 数 就 能 够 推算 出 发 、 落 点 等 预警 信 
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FE, 但 是 如 果 导 弹 的 发 点 或 落 点 位 于 高 山 、 高 原 等 
地 表 与 参考 椭 球 面 距离 较 大 的 地 区 , 使 用 该 方法 会 
降低 预警 信息 的 精度 , 且 高 程 值 越 大 , 精度 降低 得 
显然 也 就 越 多 . 因此 有 必要 改进 文献 [6] 介 绍 的 预警 
信息 推算 方法 , 使 新 的 方法 在 保持 较 高 计算 效率 和 
稳健 、 稳 定 度 的 前 提 下 尽量 消除 高 程 因 素 引 起 的 
误差 , 提高 预警 信息 的 精度 . 本 文 将 针对 这 一 问题 
展开 理论 分 析 和 探讨 . 

用 于 文 持 发 、 落 点 推算 的 高 程 信息 库 可 使 
用 全 球 数字 高 程 模 型 (Digital Elevation Model, 
DEM), 它 是 对 地 球 表面 地 形 地 貌 的 一 种 离散 数学 
表达 . DEM 的 原理 是 将 一 片区 域 划分 为 5 行 c 列 (5b、 
c 均 为 正 整 数 ) 的 四 边 形 , 计算 每 个 四 边 形 的 平均 高 
FE, 然后 以 三 维和 矩阵 的 方式 存储 这 些 平均 高 程 . 
函数 的 形式 可 描述 为 : 


V; = Gy Hi), j= 1,2,.., 


XB, Lj Wye KAA EAA, Hy AAR 
坐标 (Lj;, w;) Rb AY fe Be. OB J BN) e 2T ER AE 
球 DEM 的 生成 主要 得 益 于 全 球 卫 星 遥 感 、 卫 星 测 
高 、 船 载 测 深 等 资料 的 获取 后 . 一 般 来 说 DEME 
小 四 边 形 的 边 长 就 是 DEM 的 分 辩 率 . 常见 的 分 辨 
率 有 90 m、30 m、12.5 m. 5 m^. 分 辨 率 是 描述 
地 形 精确 程度 的 一 个 重要 指标 , 边 长 越 小 同一 片区 
域 分 割 的 四 边 形 越 多 , 这 就 意味 着 分 辨 率 越 高 , 描 
述 地 形 越 精确 . 


2 ”弹道 导弹 预警 的 原始 方法 
为 了 后 续 阐述 方便 , 简要 介绍 一 下 文献 [6] 提 出 


的 方法 . 当 导 弹 处 于 发 点 或 落 点 时 , 其 所 处 的 弹道 
位 置 与 地 球 表 面相 交 , 导弹 地 心 距 r (了 ) 与 其 星 下 点 
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A. 直接 求解 (1) 式 难度 较 高 , 因此 使 用 两 个 步骤 迭 
代 求 解 : 第 一 步 迭 代 是 在 二 体 运动 模型 下 求解 初 
始 值 ,迭代 初 值 是 由 初始 根 数 求 得 的 目标 真 近 点 
FA fo, 收敛 需要 满足 的 判定 门限 为 | 六 — Ri| < e, 其 
Pjr, — RE SS UG TUS, 目标 地 心 距 与 星 下 点 地 
心 距 之 差 的 绝对 值 , © > 0 是 收敛 门限 , 表示 第 k 次 
IER. 二 体 运 动 模型 下 , > 可 表达 为 : 
raf (1 — e°) 
1+ ecosf’ 
Ha, esp ll Fe E DIR Be B as AER 48 Fl à 
心率 , 而 有 R 可 表达 为 : 


R= s 
V1— (28 — 82) [1 — sin? isin? (w + f)] 


其 中 i、w 分 别 为 开 普 勒 根 数 中 的 轨道 倾角 和 近 地 
点 幅 角 , R, 为 地 球 极 半径 , 为 地 球 扁 率 
之 后 第 二 步 是 引入 一 阶 长 期 项 和 一 阶 短 周期 
项 的 影响 ， 对 第 一 步 处 理 获取 的 初 值 进行 修正 ， 
得 到 精确 的 交点 时 刻 . 这 一 步 又 的 收敛 门限 为 
Ir, — Ri] < e, FPR RANA UIE. XX 72035 
代 收 敛 后 可 得 到 精确 的 交点 时 刻 T0, 通过 交点 时 刻 
可 进一步 求 得 目标 在 地 固 坐 标 系 的 位 置 矢量 以 及 
轨道 与 地 球 表面 精确 交点 的 地 理 经 纬度 L 和 wy. 
3 ”有 DEM 支 持 的 精确 交点 改进 算法 
原 方法 得 出 轨道 与 参考 椭 球 体 交 点 信息 后 , 为 
了 进一步 求 得 DEM 支 持 下 更 精确 的 交点 ， 本 文 提 
出 在 原 方法 计算 结果 的 基础 上 进行 新 一 轮 的 迭代 
计算 的 新 方法 . 
在 新 一 轮 的 迭代 计算 中 ， 收 敛 的 条 件 设 置 为 
|Hp — Hc| < e2, 其 中 Hc 为 计算 求 得 的 t 时 刻 目 标 


; (3) 


S (ERRTI 
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的 高 程 , Hb 为 通过 DEM 查 询 t 时 刻 目 标 位 置 的 经 纬 
度 得 出 的 高 程 , sy > 0 是 给 定 的 门限 , 与 原 方法 的 门 
限 是 不 同 的 值 . 本 轮 迭 代 是 以 高 程 差 小 于 特定 门限 
作为 收敛 标准 , 因此 该 新 一 轮 迭 代 过 程 可 称 为 “高 
TUER WA. 具体 迭代 方法 如 下 : 

(1) 使 用 原 方法 计算 得 出 的 交点 地 理 经 纬度 L 
和 查询 DEM, SH FED. 再 通过 地 固 坐 标 系 至 
站 心地 平 坐 标 系 的 旋转 平移 转换 ,可 求 得 目标 t 时 
刻 在 以 目标 自身 为 原点 的 地 平 坐标 系 下 垂直 方向 
的 速度 vs, 求 得 时 间 修 正 量 At 为 : 


At 一 Hp /vhs 


(2) 9T =t+ At, 加 入 一 阶 长 期 项 和 一 阶 短 周 
期 项 的 影响 , 重新 外 推 计算 得 到 t 时 刻 目 标的 瞬时 
Ho, 并 由 该 根 数 计算 时刻 的 目标 地 心 距 r 和 目 
标 在 轨道 坐标 系 中 的 位 置 矢量 # 和 和 速度 矢量 六 

(3) 通 过 轨道 坐标 系 至 地 固 坐 标 系 的 坐标 旋转 
变换 可 得 到 t 时 刻 目 标 在 地 固 坐 标 系 中 的 位 置 矢量 
宛 ， 由 元 按 以 下 文献 [98] 计算 目标 轨道 与 地 球 表面 精 


lm 
s 


确 交 点 的 地 理 经 纬度 KL、 和 高 程 Hc: 
人 z(1— 6) + (26 — 6?) R,sin?o 
(1 — 6)[\/a? + y? — (26 — 62) RR. cosa ’ 
sin L = MN EN cos L — 


x 
N cosy’ N cosy’ 
Ho = Ncos^y + zsinw — Re /N, 


其 中 RR 是 地 球 赤道 半径 , a. y. 


z 分 别 为 目标 在 地 


固 坐 标 系 位 置 矢 量 的 三 个 分 量 , 此 外 : 
tana = = |a ô) 4 DURER ; 
"na r 
lal « 5. 
N= He 
vA — (26 — 82) sin? 
(48 HL. (f DEM EH Hp, JSt 
AB ps E FRR ASB nt NY ZIA UJ. E] 5 Z3 Je. HS A 
标 系 下 垂直 方向 的 速度 ws. 由 下 式 求 t 时 刻 的 修正 


*http://www.gscloud.cn/. 
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At = (Hp = Hc) [Unz * 


(5) 令 T = t+ At 回 到 本 迭代 方法 的 第 二 步 

续 计 算 tfHAt,， 回 到 本 步 又 时 ， 如 果 满 足 迭 代 条 
ed E en 
然 回 到 第 二 步 继 续 计 算 . 

迭代 收敛 以 后 , 得 到 的 计算 结果 就 是 目标 位 于 
弹道 和 地 表 的 交点 时 更 为 精确 的 时 刻 t 和 发 、 落 点 
经 纬度 LL 和 四 以 及 根据 该 经 纬度 查询 DEM 得 到 的 
FUERTE Hp. 
正如 本 文 开 头 所 述 DEM 数 据 库 中 储存 的 高 
程 点 是 离散 的 ， 这 就 意味 着 目标 轨道 与 地 表 交 点 
区 域 的 高 程 越 高 、 坡 度 越 大 、DEM 自 身分 辩 率 
越 高 时 ， 可 能 会 导致 连续 两 次 查询 得 到 的 高 程 数 
HADEN RBA, MM EIR TRE. A 
JY DTI APL, A GT DAS i eT 
限 es 的 值 ; 其 次 可 以 使 用 优选 策略 : EAE RATER 
KB kimax, “IER BOR Bl kena TRAY FE IEIE AR, XE 
TE UIER P| Ap 一 Hehi EI BS IS VCI OS 
代 计 算 求 得 的 五 、% 以 及 相应 查询 DEM 获 得 的 高 
程 Hb 作 为 最 优 近 似 值 , 视 为 最 终 计算 结果 . 


4 数值 实验 和 分 析 

由 于 缺乏 实测 弹道 数据 ,我 们 进行 了 仿真 测 
试 和 分 析 来 评估 该 方法 的 计算 效能 . 测试 中 使 用 
的 DEM 数 据 来 源 于 中 科 院 网 信 中 心地 理 数据 云 平 
台 的 SRTMDEM (Shuttle Radar Topography Mis- 
sion DEM， 航 天 飞机 雷达 地 形 测量 任务 DEM 数 
HE) 90 m 分 辩 率 高 程 数据 !， 该 数据 是 在 WGS84 
(World Geodetic System 84)itb [i] AB MADR A 
上 的 投影 . 首先 利用 己 有 的 成 熟 弹 道 设 计 软 件 设 
计 得 到 了 5 组 弹道 , 再 用 试探 法 通过 高 精度 数值 外 
推 得 到 了 这 5 组 弹道 与 地 球 表面 的 精确 交点 作为 
发 、 落 点 的 真实 值 ， 其 轨道 根 数 、 射 程 及 交点 的 
高 程 列 于 表 1. 本 文中 的 射程 是 指 惯性 空间 中 ， 导 
弹 的 轨道 在 地 球 外 部 那 一 部 分 的 弧 长 . 使 用 的 根 数 


64 卷 天 ox 学 报 5 期 


格式 是 卫星 定 轨 中 经 常 采 用 的 第 一 类 无 奇 点 根 数 : 了 仿真 计算 . 参考 文献 [6] 的 方法 , 首先 根据 弹道 根 
除 a、i 和 9 外 , 另外 的 三 个 根 数 为 5 = ecosw、7 = 数 及 发 、 落 点 选择 合适 的 观测 站 址 , 其 次 生成 使 用 
—esinu. À = M +w, 其 中 M 为 轨道 平 近 点 角 . 该 观测 站 来 跟踪 目标 的 仿真 观测 数据 , 然后 在 这 些 


在 实际 应 用 中 , 目标 身份 鉴别 及 预警 信息 生成 仿真 观测 数据 中 加 入 角度 0.1?、 距 离 50 m 的 随机 误 
是 雷达 在 跟踪 目标 时 进行 实时 处 理 的 一 部 分 , 因此 ” 25. 之 后 使 用 实验 室 已 有 的 程序 将 前 30 s 的 数据 进 
模拟 雷达 实时 处 理 过 程 , 从 中 考察 该 方法 的 计算 效 行 滤波 初始 化 , 最 后 使 用 卡尔 曼 滤 波 程序 对 30 s 后 
率 和 计算 精度 的 变化 情况 是 最 好 的 选择 . 我 们 设计 的 每 秒 一 点 的 仿真 数据 进行 滤波 来 更 新 目标 的 根 
实验 对 表 1 的 这 5 组 射程 和 高 程 各 异 的 根 数 均 进 行 HL. 


表 1 仿真 实验 中 使 用 的 轨道 根 数 、 射 程 及 发 、 落 点 高 程 数 据 


Table 1 Orbital elements, range and elevation of launching and landing points used in simulation 


experiments 
Example Time! a/km i/? Q/? € n A/?^  Range/km Height?/m Height? /m 
1 59639.8451 5544.802 159.124 24.523 0.225  —0.335 236.146 11467 4280 3699 
2 59639.8427 4999.656 72.262 198.201 —0.517 0.035 3.989 7506 2548 1488 
3 59639.8432 5158.098 32.667 285.793 0.453 0.128 164.213 7818 1132 1113 
4 59639.8404 4422.964 70.044 107.308 —0.620 —0.220 340.329 4900 370 177 
5 59639.8420 4817.056 69.139 117.594 —0.552 —0.102 349.469 6617 75 —38 


' Beijing time in the modified julian date format 
? The height of launch point 
3 The height of landing point 


每 次 滤波 更 新 根 数 后 , 我 们 分 别 采用 原 方法 和 上 的 变化 . 事实 上 , DEM 的 读 取 速 度 也 是 影响 计算 
本 文 提出 的 新 方法 产生 目标 两 个 交点 的 位 置 和 时 效率 的 一 个 重要 因素 , 但 其 涉及 的 是 数据 库 读 取 的 
Al, 观察 两 种 方法 计算 出 的 两 个 交点 时 刻 和 位 置 的 相关 技术 ,已 超出 本 文 的 讨论 范围 , 因此 仅 假 定 读 
计算 误差 变化 . 其 中 原 方法 的 收敛 精度 设 为 0.01 m, 取 DEM 的 速度 与 读 取 内 存 变量 速度 一 致 ， 不 做 进 
新 方法 的 高 程 迭代 收敛 精度 =s? 设 为 1 m, 最 大 高 程 PRF. 


— 


迭代 次 数 kwax 设 为 5 次 ,交点 时 刻 位 置 误差 的 计算 表 2 列 出 了 各 算 例 在 生成 预警 信息 时 增加 的 高 
方法 为 : 程 迭代 次 数 的 情况 , 由 表 2 可 以 看 出 ， 实验 的 5 个 算 
Ame csse TE E T BIER, VERL B MOET T TOATE BERT 

S. SPOR AUDI F REDERE, 仅 算 


H PERD EARMA kE, RAR, 例 1 和 算 例 2 的 发 点 推算 时 ， 超 过 2 次 高 程 迭 代 的 预 
是 地 固 坐 标 系 下 的 位 置 矢量 的 计算 值 和 真实 值 . 警 信息 生成 次 数 较 多 , 甚至 未 收敛 的 生成 次 数 分 别 
至 于 计算 效率 方面 的 变化 , 新 方法 的 计算 量变 为 72 次 与 99 次 , 这 是 由 于 这 两 个 发 点 处 于 高 程 高 、 
化 主要 来 源 于 高 程 迭 代 所 增加 的 迭代 计算 , 因此 在 坡度 大 的 山地 中 . 算 例 1 的 落 点 位 于 玻利维亚 的 乌 
新 方法 计算 效率 的 评价 方面 , 主要 观察 增加 的 高 程 尤 尼 盐 沼 , 此 处 高 程 变化 很 小 , 因此 虽然 该 点 高 程 
迭代 次 数 的 多 少 . 假定 滤波 过 程 中 有 i 次 预警 信息 达 3799 m, 却 仍然 能 保证 大 体 上 和 迭代 2 次 即 收敛 . 至 
生成 时 进行 了 m (m = 1,2,…:) 次 高 程 迭 代 , mi 于 发 、 落 点 高 程 都 很 低 的 算 例 5， 全 程 高 程 迭代 基 
值 可 以 用 来 评价 新 方法 相 比 于 原 方法 在 计算 效率 本 1 次 收敛 . 


=| 
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表 2 新 方法 生成 预警 信息 时 增加 的 高 程 迭代 次 数 的 情况 


Table 2 Increase of elevation iterations when the new method generates warning information 


Launch Point Landing Point 
Example 
i? 2 3 4 5 NC? 1 2 3 4 B NC 
1 0? 12 43 94 49 72 13 235 20 1 0 1 
2 0 108 57 6 0 99 0 213 48 2 0 7 
3 48 222 0 0 0 0 5 244 18 0 0 3 
4 0 270 0 0 0 0 0 265 3 0 0 2 
5 225 44 0 0 0 1 270 0 0 0 0 0 


! The number of elevation iterations increased by the new method (i.e. the l in the text) 
? Unconvergence 
3 This value refers to the number of times of generating early warning information for the elevation iteration 


with the number of times shown at the head of the column increased (i.e. the m in the text) 


图 1、3、5、7、9 展 示 了 5 个 算 例 的 弹道 目标 参考 椭 球 的 算 例 5,， 两 种 方法 推算 的 交点 位 置 和 交 
发 、 落 点 位 置 误差 随 雷达 跟踪 时 间 的 变化 情况 ， 点 时 刻 差 距 很 小 , 对 于 高 程 75 m 的 发 点 , 新 方法 精 
图 2、4、6、8、10 展 示 了 5 个 算 例 弹 道 目 标 发 、 落 。 ” 废 仅 略微 优 于 原 方法 , 而 对 于 高 程 仅 -38 m 的 落 点 ， 
点 时 刻 误 差 随 雷达 跟踪 时 间 的 变化 情况 . 由 图 1 至 图 9 和 图 10 显 示 两 种 方法 的 计算 结果 曲线 重合 在 
图 10 可 以 看 出 , 在 滤波 初始 阶段 ,由 于 定 轨 误差 较 。” 起, 意味 着 两 种 方法 的 推算 精度 几乎 一 致 


大 ,因此 新 方法 和 原 方法 计算 的 时 刻 和 位 置 误差 以 上 计算 结果 表明 , 交点 位 置 高 程 值 越 大, 使 

D c FETE AY LIES. 交点 时 刻 和 TAE i A 
AMER LENE, AMEERDATEAE ARMS RUIT RR UG AM 
KG, ANAE ERATES OEMS. EA Ne, 
ASHES WENERTT3 54 mung, ean 全 
刻 和 位 置 误差 均 保持 稳定 , 此 时 可 以 看 到 , 对 于 前 AMOS, SUERTIE EN ERAMI- 
ORARE TARARE 新 方法 计算 得 由 。 2 次 高 程 迁 代 的 计算 量 在 交点 高 程 高 、 坡 度 大 时 | 
三。 的 交点 时 刻 和 交点 位 置 精度 都 好 于 原 方法 . 交点 位 迁 代 次 数 会 增加 , 未 收敛 的 情况 会 增多 , BE 


C 0 显 .对 于 算 例 ! 的 弹道 交点 时 刻 误差 ， 新 方法 的 精 。 ”入 高 程 迁 代 后 ,可 在 增加 少量 迁 代 次 数 的 情况 下 ， 
度 为 0.5 s 左 右 , 交点 位 置 误差 为 3.8 km 左右 , 而 原 极 大 程度 地 消除 高 程 带 来 的 影响 . 相 比 于 其 他 算 例 ， 
方法 给 出 的 两 种 误差 分 别 为 2 s 左 右 和 10 km 左右 ， ” 算 例 1 采用 新 方法 的 计算 精度 也 偏 低 . 这 是 因为 其 
对 于 算 例 2-4, 新 方法 的 交点 位 置 和 交点 时 刻 误差 。 主要 误差 来 源 是 较 长 射程 导致 的 更 久 的 外 推 时 间 ， 
分 别 为 百 米 量 级 和 0.01 s 左 右 ， 均 优 于 原 方法 公里 进而 使 得 测 轨 误 差 有 更 长 时 间 的 放大 . 而 这 一 部 分 
量 级 和 0.02 s 左 右 的 误差 . 对 于 发 、 落 点 位 置 接近 ”的 计算 误差 显然 是 无 法 通过 高 程 迭 代 消 除 的 . 
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图 1 两 种 方法 对 算 例 1 的 发 、 落 点 位 置 计算 误差 随 跟踪 时 间 的 变化 情 “图 2 两 种 方法 对 算 例 1 的 发 、 落 点 时 刻 计算 误差 随 跟踪 时 间 的 变化 情 
况 Ui 
Fig.1 'The variation of the calculation error with the Fig.2 The variation of the calculation error with the 
tracking time for the launch and landing distance of example tracking time for the launch and landing time of example 1 
1 by the two methods by the two methods 
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图 3 两 种 方法 对 算 例 2 的 发 、 落 点 位 置 计 算 误差 随 跟踪 时 间 的 变化 情 “图 4 两 种 方法 对 算 例 2 的 发 、 落 点 时 刻 计算 误差 随 跟 踪 时 间 的 变化 情 
D 况 
Fig.3 The variation of the calculation error with the Fig.4 The variation of the calculation error with the 
tracking time for the launch and landing distance of example tracking time for the launch and landing time of example 2 
2 by the two methods by the two methods 
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Fig.5 The variation of the calculation error with the 
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算 例 3 的 发 、 落 点 位 置 计算 误差 随 跟 踪 时 间 的 变化 情 图 6 
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Fig. 7 The variation of the calculation error with the 


tracking time for the launch and landing distance of example 
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图 7 两 种 方法 对 算 例 4 的 发 、 落 点 位 置 计 算 误差 随 跟踪 时 间 的 变化 情 图 8 
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Fig.6 The variation of the calculation error with the 
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Fig.8 The variation of the calculation error with the 
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tracking time for the launch and landing time of example 3 


两 种 方法 对 算 例 4 的 发 、 落 点 时 刻 计算 误差 随 跟踪 时 间 的 变化 情 


tracking time for the launch and landing time of example 4 
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图 9 两 种 方法 对 算 例 5 的 发 、 
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Fig.9 The variation of the calculation error with the 


tracking time for the launch and landing distance of example 
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5 ”结论 和 讨论 


本 文 基于 文献 [6] 的 研究 提出 了 一 种 有 高 程 信 
息 库 DEM 支 持 的 预警 信息 生成 方法 , 仿真 实验 结 
果 表 明 , 该 方法 具有 下 列 特点 : 


(了 该 方法 相 较 于 原 方 法 仅 增 加 了 少量 高 程 迭 
代 计 算 ， 且 对 个 能 收敛 的 情况 有 截断 和 优选 措施 ， 
因此 能 保持 较 高 的 计算 速度 、 稳 健 和 稳定 度 , 能 够 
Vol o Sac 警 信息 , 以 供 导弹 防御 决策 

(2) 在 雷达 跟踪 弹道 导弹 的 实时 处 理 过 程 中 , 该 
方法 计算 的 预警 信息 相 比 原 方法 能 够 提升 交点 位 
置 和 交点 时 刻 的 预报 精度 . 交点 位 置 的 高 程 值 越 大 ， 
情 度 提升 越 明显 . 通常 在 滤波 进行 了 3.5-4 min 后 ， 
能 够 提供 比较 准确 的 发 、 落 点 信息 ， 发 、 落 点 位 
置 和 时 刻 的 提升 幅度 可 分 别 达 到 6 km 和 1.5 s. 对 于 
发 、 落 点 高 程 值 越 小 的 弹道 目标 , 其 发 、 落 点 的 计 
算 精 度 的 提升 效果 会 降低 . 

此 外 , 本 文 设 计 的 仿真 实验 显示 , 在 设备 观测 
的 前 2-3 min, 原 方法 和 新 方法 计算 的 精度 都 不 高 ， 
甚至 出 现 了 新 方法 计算 的 预警 信息 精度 低 于 原 方 
法 的 情况 . 因此 在 实际 应 用 中 , 可 以 考虑 在 滤波 达 
到 稳 态 之 前 仅 使 用 原 方法 进行 目标 身份 鉴别 和 生 
成 预警 信息 , 待 滤波 达到 稳 态 后 再 加 入 高 程 迭代 计 


落 点 位 置 计算 误差 随 跟踪 时 间 的 变化 情 
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Fig.10 The variation of the calculation error with the 


tracking time for the launch and landing time of example 5 


by the two methods 


算 , 这 样 就 可 以 在 不 影响 总 体 预警 信息 精度 的 前 提 
F, 减少 多 次 高 程 迭 代 计 算 以 及 DEM 数 据 库 读 取 
操作 , 从 而 在 观测 设备 跟踪 过 程 中 进一步 提高 计算 
速度 和 效率 . 

需要 说 明 的 是 , 在 实际 情况 下 目标 在 主动 段 的 
受 力 过 程 难以 获知 , 因此 即使 使 用 本 文 的 方法 计算 
的 发 点 信息 仍然 会 有 较 大 误差 . 落 点 计算 则 不 同 ， 
设计 弹道 一 般 要 求 在 不 考虑 大 气动 力 影响 的 情况 
下 准确 到 达 落 点 . 在 此 阶段 的 控制 调整 主要 任务 是 
使 目标 能 准确 击 中 指定 地 点 , 因此 , 采用 本 文 的 方 
法 来 提升 落 点 位 置 时 刻 的 精度 会 有 较 大 的 应 用 价 
值 . 


参考 文献 


TEE. 弹道 导弹 中 段 雷达 
天 大 学 , 2020 

AÈ, 何 广 军 , 高 鹏 . 信息 通信 , 2017, 9: 113 

RE. 弹道 导弹 中 段 目标 识别 研究 . 南京 : 南京 大 学 , 2012 

魏 五 洲 , EF, PAWS. 兵 工 自动 化 , 2022, 41: 70 

KER, TER, 黄金 才 等 . 弹道 学 报 , 2012, 24: 22 

Fa, RI, 陈 务 深 . 天 文学 报 , 2014, 55: 256 

FIRR, HS, TKF. 武汉 大 学 学 报 (信息 科学 版 )，2018，43: 
1927 

8] Bowring B R. Survey Review, 1985, 28: 202 


um 


标识 别 方法 研究 . 南京 : 南京 航空 航 


s åO a ROC LM 


64 卷 杨 冬 等 : 有 高 程 数据 库 支 持 的 弹道 预警 信息 生成 方法 5 期 


The Generation of Trajectory Early Warning Information with 
Elevation Database 


YANG Dong? XU Jint? DU Jian-li!? 


(1 Purple Mountain Observatory, Chinese Academy of Sciences, Nanjing 210023) 
(2 Key Laboratory of Space Object and Debris Observation, Chinese Academy of Sciences, Nanjing 210023) 


ABsTRACT This paper proposes a new method of deriving the launch and landing locations based on the 
orbital motion characteristics of the missiles and elevation database. In this method, the DEM support is 
introduced to change the assumption of the earth from the reference ellipsoid to a more real shape. The 
method to improve the prediction accuracy of the position and time of the launch and landing points is 
discussed. The results with DEM (Digital Elevation Model) indicate the accuracy can be improved by 
adding a few elevation iterations on the original method. The higher the elevation of the starting and 
landing points, the better the result would be. 


Key words ballistic trajectory, methods: numerical, height 


